Manifestation of the Roughness-Square-Gradient Scattering 
in Surface-Corrugated Waveguides 
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We study a new mechanism of wave/electron scattering in multi-mode surface-corrugated waveg- 
uides/wires. This mechanism is due to specific square-gradient terms in an effective Hamiltonian 
describing the surface scattering, that were neglected in all previous studies. With a careful analysis 
of the role of roughness slopes in a surface profile, we show that these terms strongly contribute to 
the expression for the inverse attenuation length (mean free path), provided the correlation length 
of corrugations is relatively small. The analytical results are illustrated by numerical data. 

PACS numbers: 42.25.Dd; 73.21.Hb; 73.23.-b;73.50.Bk;73.63.Nm 



I. INTRODUCTION 

The subject of wave transport through guiding surface- 
disordered systems is of great importance in various phys- 
ical applications. Recently, this topic has attracted even 
more attention due to a burst of developments in nano- 
science where frequently one deals with devices or struc- 
tures (quantum wires, leads, etc.) whose surface's irreg- 
ularities become more important than those in the bulk. 
Other situations for either electromagnetic or acoustic 
waves (remote sensing, photonic and acoustic devices, 
optical thin films, etc.) are analogous. Therefore, in 
what follows we do not make any distinction between 
electromagnetic/electron or any other type of waves. 

In spite of an extensive research of unitary wave scat- 
tering from rough surfaces, the problem of transport in 
the waveguides with such profiles still remains open. This 
problem is a great theoretical challenge since it deals with 
the multiple scattering of a wave from lateral walls. As 
a result of this scattering, the unperturbed longitudinal 
wave number k n of an n-th propagating normal mode 
changes, k n — > k n + Sk n , due to a complex amount 5k n , 



the theory accounting for the guided and scattered waves, 
is not a trivial task. Many approaches have been pro- 
posed i n co n nect ion with the surface scatteri ng ( see, e.g., 
Refs. Ill l2l3l l4l l5l6l l7ll8l l9ll0lll ill 1 2ll 1 3ll 1 4ll 1 5ll 1 6ll 1 7ll 1 SI and ref- 
erences therein). One of the main tools to treat this 
problem is a reduction of the surface scattering to the 
bulk one in such a way that the latter can be formally 
described by an effective Hamiltonian H, 



(1.2) 



Sk n = 7„ + i(2L n )~ 



(1.1) 



The real part j n is responsible for the roughness-induced 
correction to the phase velocity, while L n is the elec- 
tron total mean free path, scattering length or attenuation 
length of a given mode. As is known, the shift 7„ does 
not change the static transport of a disordered systems. 
Therefore, our further analysis shall be focused only on 
the attenuation length L n . 

Evidently, statistical properties of a rough surface pro- 
file £(x) give a strong impact on the scattering process. 
However, a proper incorporation of these properties into 



with flat boundaries, however, with a complicated po- 
tential U. To the best of our knowledge, originally the 
idea of this approach was discussed by Migdali^. Af- 
ter, it was frequently used in the theories of c lass ical 
and quantum wave scattering, see, e.g., Refs. I^T^ l8l9lll0l 
IT1II12II 13113 One should stress that in the majority of 
them2i2i2iiSiii the study was restricted to the lowest or- 
der in the root-mean-square roughness height a . Other 
methodsi2ii&i4ii5ii£iiLi£ were mainly based on a princi- 
pal assumption that the surface roughness is sufficiently 
smooth. 

In this paper we present analytical results demonstrat- 
ing a new scattering mechanism missed in previous stud- 
ies of the surface scattering. This mechanism is due to 
specific square- gradient (SG) terms that are proportional 
to a 2 £,' 2 (x) in the potential U. We argue that the discov- 
ered scattering mechanism is substantially different from 
the known ones associated with the roughness amplitude 
and roughness gradient of the surface. 

The two last mechan isms have been already studied 
(see, e.g. Refs. Hill fi 111) . Their main contributions were 
shown to depend on the terms containing the quantities 
a^(x) and a£'(x), respectively. In the performed analysis 
the square-gradient terms related with the new scatter- 
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ing mechanism were discarded due their seemingly small 
contribution. Indeed, the square-gradient terms are for- 
mally proportional to a 2 , similar to other terms that arise 
in the next (second-order) approximation in the ampli- 
tude and gradient roughness. However, we have found 
that the square-gradient terms have a very strong depen- 
dence on the roughness correlation length R c , in contrast 
with those terms appearing in the second-order approx- 
imation. Specifically, with a decrease of R c the square- 
gradient terms in the expression for the mean free path 
compete with the standard terms proportional to crt;(x) 
and <r£'(x). 

One should stress that our approach is restricted by 
the first order approximation, with a careful analysis of 
all terms that may be important in this approximation, 
and taking into account an unexpectedly strong influence 
of the square-gradient terms. Our goal is to study the 
new scattering mechanism and establish the conditions 
when it should not be neglected. For this, we derive 
a correct expression for the attenuation length L n that 
incorporates the contribution of the above mechanism. 
One should emphasize that the approach we use does not 
assume any special restrictions to the model parameters 
(such as the smoothness of surface profiles) except for 
general conditions of weak scattering. 

The paper is organized as follows. In Sec. ^ we formu- 
late the problem and discuss the coordinate transforma- 
tion used to represent the surface scattering as the bulk 
one. Our approach involves an average Green's func- 
tion whose longitudinal wave number is a modification of 
the unperturbed one. Also, all expressions corresponding 
to the unperturbed problem are given and discussed in 
this Section. In Sec. IIIII we derive a Dyson- type integral 
equation for the exact Green's function. From the ex- 
act expression for the scattering potential we design two 
Hermitian random operators. The first one is associated 
with both the amplitude (AS) and the gradient scatter- 
ing (GS) mechanism, while the second one is associated 
only with the square- gradient scattering (SGS) mecha- 
nism. The former operator gives rise to the roughness- 
height power (RHP) spectrum W(k x ), while the latter to 
the roughness-square-gradient power (RSGP) spectrum 
T[k x ). 

In Sec. El we obtain the averaged Green's function 
by applying a perturbative method with respect to the 
above operators. Since in this work we restrict ourselves 
to the analysis of the attenuation length L n of the n-th 
conducting mode, we focus our attention to the imagi- 
nary part of the proper self-energy. In Sec.[V] in corre- 
spondence with a clear independence between AS+GS 
and SGS mechanisms, we develop a natural approach 
to define two attenuation lengths, the well known one, 
and the new SGS length l\? . These two lengths 
are related to the RHP and the RSGP spectrum, and, 
as a consequence, they are characterized by a 2 and a 
dependencies, respectively. We study their interplay for 
two limit cases of the roughness surface, for small-scale 
and large-scale roughness. We discuss here the main re- 



sult according to which the larger roughness slope (t/R c , 
the larger contribution of the SGS mechanism. We also 
present a numerical analysis assuming that the surface 
profile £(x) has the standard Gaussian binary correlator. 
Finally, in Section lvTl we outline our conclusions. A short 
presentation of main results can be found in Ref. l20i 



II. PROBLEM FORMULATION 

In what follows, we consider an open plane waveg- 
uide (or conducting quasi-one-dimensional wire) of av- 
erage width d, stretched along the ir-axis. For simplicity, 
one (lower) surface of the waveguide is assumed to be 
flat, z = 0, while the other (upper) surface has a rough 
profile, z — d + <j£(x), with a as the root-mean-square 
roughness height. In other words, the waveguide occupies 
the region 

-co < x < oo, < z < w(x) (2.1) 

of the (x, z)-plane. Here the fluctuating wire width w(x) 
is defined by 

<w(x) = d + a£(x), (w(x))=d. (2.2) 

The random function £(x) describes the roughness of 
the upper boundary. It is assumed to be a statisti- 
cally homogeneous and isotropic Gaussian random pro- 
cess with the zero mean and unit variance, 

(£(*)) =0, (£ 2 (x))=l, {S(x)t(3>))=W(x-x>). 

(2.3) 

Here the angular brackets stand for statistical averaging 
over different realizations of the surface profile £ (x) . We 
also assume that its binary correlator W(x) decreases on 
the scale R c and is normalized to one, W(0) = 1. 

The roughness-height power (RHP) spectrum W(k x ) 
is defined by 

/oo 
dxexp(-ik x x)W(x), (2.4) 
-oo 

Since W(a;) is an even function of x, its Fourier transform 
(|2.4|l is even, real and non-negative function of k x . The 
RHP spectrum has maximum at k x — with W(Q) ~ R c , 
and decreases on the scale R~ l . 

In order to analyze the surface scattering problem we 
shall employ the method of retarded Green's function 
G(x, x'; z, z'). Specifically, we start with the Dirichlct 
boundary-value problem 

( d 2 d 2 \ 

{d^ + d^ 2+k J K x >**>*) 

= S(x-x')S(z- z'), (2.5a) 

G(x, x ; z = 0, z ) = 0, 

g(x,x';z = d + a£{x),z') = 0. (2.5b) 
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Here S(x) and S(z) are the Dirac delta- functions. The 
wave number k is equal to uj/c for an electromagnetic 
wave of the frequency w and TE polarization, propagat- 
ing through a waveguide with perfectly conducting walls. 
For an electron quantum wire, k is the Fermi wave num- 
ber within the isotropic Fermi-liquid model. 



A. Coordinate transformation 

The equation l|2.5a|) does not contain any scattering 
potential. In contrast with the bulk scattering, the elec- 
tromagnetic/electron waves experience a perturbation 
due to scattering from the upper wall, therefore, the per- 
turbation is hidden in the boundary condition <|2.5b[) . In 
order to formally describe the surface scattering as a bulk 
one, we perform the canonical transformation to new co- 
ordinates, 



Xold i 

d 

w(x) 



Zold = 



d + er£(x 



■Zold, 



(2.6a) 
(2.6b) 



in which both waveguide surfaces are flat. Correspond- 
ingly, we introduce the canonically conjugate Green's 
function (below for convenience we drop the subscript 
"new" for x and z), 



. . v/ w(x)w(x') . . A , . 

Qnew{x, x;z,z) = Qoidix, x ;z,z). (2.7) 

a 

As a result, we arrive at the equivalent boundary- value 
problem governed by the equation 



d 2 



dx 2 dz 2 



d 2 \ , 
+ -5-0 + k 2 G(x,x';z, z) 



-U{x, z) G(x, x'; z, z') = 6(x - x')5(z - z'), (2.8a) 



G(x,x';z = 0,z') = 0, 
Q(x, x'; z = d, z ) = 0. 



(2.8b) 



B. Unperturbed Green's function 

The unperturbed Green's function (?q(|x — x'\;z,z') 
obeys the boundary- value problem l|2.8[l with U(x, z) = 
(when cr = 0, and therefore, w(x) = d). It is determined 
as follows, 



Go(\x - x'\;z,z') = 



exp \i[k x x — k x x )\ 



(27T) 



x 90 (k x X)- n{kzZ)MKZ,) 



k z k' z 



(2.10a) 



N d 

sm(— jsm 



imz\ . ( nnz \ exp (ik n \x — x'\) 



(2.10b) 



Here k x and k z are the lengthwise and transverse wave 
numbers, 



k z = k z (k x ) = Vk^k 2 , k' z = k z (k' x ). (2.11) 
Their unperturbed eigenvalues are k n , 



k n = Vk 2 - (nn/d) 2 , n = 1, 2, 3, . . . , N d , (2.12) 

and im/d respectively. The total number Nd of propa- 
gating waveguide modes (conducting channels that have 
real values of k n ) is determined by the integer part [. . .] 
of the ratio kd/n, 



N d = [kd/ir] 



(2.13) 



The first expression (|2. l()a,|) gives the unperturbed 
Green's function in the pole representation where 
go(k x , k' x ) is called the pole factor, 



9o(k x , K) = 2ir8(k x - k x ) g (k x ), 
9o(k x ) = k z cot(k z d) -> 



N d 

E 



(irn/d) 2 



kfid V kx ~\~ kfi ~\~ ^0 kx kn zO 



(2.14a) 
(2.14b) 

. (2.14c) 



Here the effective surface scattering potential U(x, z) is 
given by the following expression, 



U(x,z) = 



1 - 



w 2 (x) 



dz 2 



w(x) 

a 2 e 2 {x) 
w 2 {x) 



C'(x) + 2^'(x)- 



i d 

2 +Z Yz 



-+3z— + z 
4 dz 



dz 



(2.9) 



Note that the prime over the function £(x) denotes a 
derivative with respect to x. It is important to stress 
that Eqs. I|2.8|l - (|2.9|) are exact and valid for any form of 
the surface profile £(x). 



The evaluation of the integrals in Eq. (|2.10a|) over the 
poles of go(k x ,k' x ), results in the mode representation 
(j2.10b|) for the unperturbed Green's function. 



III. DYSON EQUATION 

With the use of the Green's theorem, it can be easily 
shown that the boundary-value problem (|2.8|l is directly 
reduced to the following exact Dyson-type integral equa- 
tion 

/oo I'd 
dx\ I dzi 
-oo J0 

G (\x - a;i|; z, zi)U(xi,zi) Q(x 1 ,x'\ z\, z'). (3.1) 
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This equation relates the perturbed by surface disor- 
der Green's function Q(x, x'\ z, z') to the Green's func- 
tion Q (\x — x'\; z, z') of the waveguide with perfectly flat 
boundaries. 



Similarly to the pole representation (|2.10a(l for the un- 
perturbed Green's function, let us seek Q(x, x') z, z') in 
the form 



Q(x,x';z,z') = 



rj- exp [i(k x x - k x x )\ 



(2tt)2 



wz V \ sin(fc 2 z)sin(/c^') 
Xg{k x ,k x ) - 



k z k' z 



(3.2) 



for the perturbation potential, 



a(k x , k x ) 



w 2 (xi) 



dxi f d dzi expH^)^ 1 

JO *z 

dz\ 



+i(k x+ k' x )^. 



1 d 



<tV 2 (si) 
w 2 {x\) 



x exp (ik' x xx) 



4 dzi dz{ 

sa\(k' z z\) 



(3.5) 



With this method the problem of deriving the Green's 
function Q(x, x'\ z, z') is reduced to obtain its pole fac- 
tor g(k x ,k' x ). To this end, we substitute the pole rep- 
resentations <|2.1Ua|) and (|3.2|) into Eq. (|3.1|l and get the 
non-local Dyson equation in the fc^-representation, 



g(k x ,k' x ) = g (k x ,k' x ) 
dq x dq' x 



{2nf 



)S(q x ,q x )g(q x ,k x ). 



(3.3) 



In this representation the effective surface scattering po- 
tential S(fe x , k' x ) is 



_ , , , / / , , sin(fc z zi) 
r.{k x ,k x )= I dxi I dzi exp {-ik x xi) 

sm(k z zi) 



xU(xi,z{) exp(ik' x xi) 



(3.4) 



Thus, we obtained the Dyson- type integral equa- 
tion 13.3fl with an exact expression l|3.4|) for the kernel 
E,(k x ,k' x ). After the substitution of the expression (|2.9|) 
for U, we realize that the kernel consists of three groups 
of terms. The first one has the factor (1 — d 2 /w 2 (xi)) 
while the second and third group has, respectively, the 
factor a/w{x\) and a 2 jw 2 {x\). One needs to note that 
while the kernel S(k x ,k' x ) is Hcrmitian in the whole, its 
latter two parts are non-Hermitian individually. To have 
each part of E,(k x ,k' x ) Hermitian, we perform the inte- 
gration by parts for the term containing £"(xi) in the 
second group. After that we come to the final exact form 



This equation has a peculiar structure very useful in fur- 
ther analysis. The kernel written in this form also con- 
sists of three groups of terms, however, now they repre- 
sent different scattering mechanisms. 

Since we are interested in the averaged Green's func- 
tion, we have to calculate the binary correlator of 
E(k x ,k' x ). Therefore, to avoid very cumbersome calcula- 
tions, it is reasonable to make a simplification of , k' x ) 
that does not destroy the Hermitian structure of each 
group of terms. To this end, and taking into account 
that our interest is in the typical case of small surface 
corrugations (a <C d), we can do the following: expand 
the factor 1 — d 2 /w 2 (xi) « 2a^{x\)/d in the first term 
of Eq. (|3.5|) and put w{x\) rs d in all the others. In such 
a way, we get a suitable approximate expression for the 
surface scattering potential, 



/oo r-d 
dxi I dzi exp(-ik x Xi) 
-oo J0 



sa\{k z zi) 



2o£(xi) d 2 



d dzi 



+ i(k x + k' x ) 



d 2 



x exp (ik' x xi) 



d 

1 „ 2 d 2 
- + 2z 1 — + z 2 — 1[ 
4 dzi oz{ 



1 d 

2 +Zl dz~ 1 



sin(k' z zi) 



(3.6) 



The expression H3.6(l contains one term that depends 
on the amplitude of the roughness profile a^(x±), and 
two groups of terms that depend on the roughness gradi- 
ent o~£'(xi) and the roughness square- gradient a 2 ^' 2 (xi) 
terms, respectively. An interesting point to mention is 
that the last group, due to its proportionality to a 2 , was 
neglected in all previous studies of transport properties 
of surface-disordered waveguides. However, as we show 
below, the scattering due to these terms has the proper- 
ties very different from those described by other terms, 
and should be properly taken into account. 
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In order to proceed further, we assert that the kernel 
can be written as the sum of its average, (E.(k x , k' x )), and 
fluctuating, 3(Je x ,k' x ), parts, i.e., 



a(fc x , k x ) — (a(k x , k x )) + ^{k x , k' x ). 



(3.7) 



It can be tested that the average, (H(fc x , k' x )), contributes 
only to the real part j n of the complex renormalization 
8k n of the lengthwise wave number k n (see Eq. I|l.l|l ). and 
therefore, does not change static transport properties of 
the surface disordered waveguide. Thus, we will omit it 
because our interest is in the attenuation length L„, not 
in 7 n . 

From Eq. 1)3. 6)1 one can see that all terms with the 
linear dependence on a have zero mean-value. There- 
fore, the average part (S(k x , k' x )) of the kernel S(fc x , k' x ) 
is associated only with the group of terms containing 
a 2 t;' 2 (xi). To extract their fluctuating contribution, wc 
should subtract from £' (xi) its mean value, (£' In 
such a way, we introduce the following zero-mean-valued 
operator 



V(*)=£' 2 (*)-(£' 2 (*)>, 



(V(x)} = 0. 



(3.8) 



The operator V{x) plays a special role in our further con- 
sideration. In accordance with the Gaussian nature of the 
surface-profile function £(x), the operator V(x) is uncor- 
related with both and £'(x), 

(Z(x)V(x'))=0, (e(x)V(x'))=0. (3.9) 
Its pair correlator is given by 

(V(z)V(x')) = 2(£'(x)£V)> 2 = 2W" 2 (x - x'). (3.10) 

Since we use the fc^-representation, it is worthwhile to 
define the Fourier transform of the operator V(x), 



V(k x ) 



dxexp(— ik x x)V(x). (3.11) 



Also, we need the correlator of its Fourier transform, 

(V(k x )V(k' x )) = An8(k x + k' x ) T(k x ). (3.12) 

Here the roughness-square-gradient power (RSGP) spec- 
trum T(k x ) is 



T{k x ) 



dxexp (— ik x x)W" 2 (x) 



(3.13) 



One should stress that on the one hand, through the 
integration by parts, the power spectrum of the rough- 
ness gradients can be reduced to the RHP spec- 
trum W(k x ). On the other hand, it is not possible to do 
the same for the RSGP spectrum T(k x ). This very fact 
reflects a highly non-trivial role of the square-gradient 
scattering (SGS), giving rise to its competition with the 
well known scattering mechanism, in spite of the seeming 
smallness of the term <7 2 £' 2 (x). 



With the introduction of V(k x ) we are ready to explic- 
itly write down the fluctuating part of the total scattering 
potential, 



E(k x ,k' x ) = Ei(k x ,k' x ) + E 2 (k x ,k' x ). 



(3.14) 



The first summand is associated with the first and second 
terms in Eq. 1)3. 6)1 . 

- „ f°° , f d , , s sin(fc z zi) 

ai(k x ,k x ) = / dxi / dzi cxp(-ifc a; xi) 

J — oo JO 



2a£,( Xl ) d 2 



i(k 



d dz\ 

sin(fc^i) 



x T- K x ) j 



1 



k z 

d 



x exp{ik x x 

i\, z 

= -a[A(k x ,k' x ) + B(k x , k' x )]l{k x - k' x ). 



(3.15a) 



(3.15b) 



The second summand is related to the third term in 
Eq. H3.6fl . with £ /2 (xi) being replaced by V(xi), 



/oo rd 
dxi I dziexp(-i/c 3 
-oo Jo 



Xx) 



sin(fc z zi) 



. t 2 V(si) 
d 2 



4 +2 *5£ 



x exp(ik x xi) 



sm(k' z z{) 



= a 2 D(k x ,k' x )V(k x -k x ). 



(3.16a) 



(3.16b) 



In Eq. (|3.15(l the quantity £(k x ) is the Fourier trans- 
form of the function £ (x) , 



dx exp(— ik x x) (3-17) 



Also, we have introduced the following quantities in the 
first summand: 



Ik! f d 

A(k x ,k' x ) = -j-^ J dz sm{k z z)sm.{k' z z), 

A(^nk n , zhfc n ') = 8 nn i . 



(3.18a) 
(3.18b) 



d_ 

dz 



sm(k' z z) 



(3.19a) 
(3.19b) 



B(±k n ,±k n t) = cos[7r(n - n')](l - 6 nn >). 
In the last expression we have used the equality 

k 2 x -k' x 2 = -{kl-k' z 2 ) (3.20) 



6 



that directly follows from the energy conservation k 2 = 
k^+k 2 = k' x 2 +k' z 2 (see the definition fFTfy ). And finally, 
in the second summand the following factors appeared, 



7 + 2z 7T + z 
4 dz 



D(±k n ,±k n ,) 



2d 



dz 2 

1 

3 



sm(k z z) 
k z 

sin(fc^z) 
k' z 

1 



(3.21a) 



7r 2 (n 2 — n' 2 )' 



(27m) 2 _ 

cos[tt(?i - n')](l - <W)< (3.21b) 



The operator E(k x ,k' x ) written as the sum l|3.14|l of 
specially designed terms, has a very convenient form. 
First, both terms are chosen to have zero average. Sec- 
ond, since the operators £(x) and V{x) do not correlate 
with each other (see Eq. (|3.9|l ^. their Fourier transforms 
are also uncorrelated, 



m,Mk'j) = o. 



(3.22) 



Due to the condition (|3.22|) . the scattering potentials 
Si(k x ,k' x ) and E^fcc, k' x ) are also uncorrelated. How- 
ever, they have the following autocorrelators 

(Zx(k x ,q x )'Zi{q x ,k' x )) = 

2TT6(k x -k' x )Q 1 (k x ,q x ), (3.23a) 

Qi(k x ,q x ) = a 2 W(k x - q x ) [A(k x ,q x ) + B(k x ,q x )] 

x[A(q x ,k x )+B(q x ,k x )]. (3.23b) 



(E 2 (k x ,q x )E 2 (q x ,k x )) = 

2nS(k x -k' x )Q 2 (k x ,q x ) 7 (3.24a) 



Q2(k x ,q x ) = 2a 4 T(k x - q x ) 
xD(k x ,q x )D(q x ,k x ) 



(3.24b) 



As a result, the correlator of the fluctuating scattering 
potential l|3.14|) in the fc^-representation, is as follows 

(E(k x ,q x )E(q x , k' x )) = 2nd(k x - k' x )Q(k x , q x ), (3.25a) 

Q(k x ,q x ) = Qi(k x , q x ) + Q2(k x ,q x ) ■ (3.25b) 

It is now possible to perform an appropriate pertur- 
bative averaging of the Dyson equation (|3.3I) , and at the 
same time, to separate a relative contribution of the SGS 
mechanism from the total scattering process. 

IV. AVERAGE GREEN'S FUNCTION 

Now we are in a position to replace the problem for the 
random Green's function G(x, x'; z, z') with the problem 



for the Green's function (G(x, x'; z, z 1 )) averaged over the 
surface disorder. It is evident that (G(x, x'; z, z')) is gov- 
erned by Eq. (|3.2|l with the average pole factor (g(k x , k' x )) 
instead of the random one. To perform the averaging of 
Eq. 13.3f> with E(k x ,k' x ) given by Eq. I|3.14|) and obtain 
{g{k x ,k' x )) , we can apply one of the standard and well 
known perturbative methods. For example, it can be the 
diagrammatic approach developed for surface disordered 
systems!, as well as the technique developed in Ref. 3. 
Both of the methods allow one to develop the consistent 
perturbative approach with respect to the scattering po- 
tential, that takes adequately into account the multiple 
scattering from the corrugated boundary. Then, we come 
to the following result 



(g(k x ,k' x )) = g a (k x ,k' x 
dq x dq' x dq"dq 



(2nY 



9o(k x ,q x ) 



x (-(Qx,q x )go{q x ,q' x h(q x ,q x '))(g(q x ',k' x )). (4.1) 

Taking into account the presence of delta-functions in 
the definitions l|2.14a(l for the unperturbed pole factor 
and in Eq. (|3.25a|) for the scattering potential, we can 
take explicitly the integrals over q x , q x and q x in the 
second term. After that, Eq. (|4.1I) becomes an algebraic 
one. Its solution has the form, 



(g(k x ,k' x )) = 2n6(k x - k' x ) g(k x ) 



(4.2a) 



9(k x ) = [g^{k x )-M{k x )\ \ (4.2b) 

The quantity M(k x ) is called the self-energy or mass op- 
erator. It is described by the formula 

-7rQ{k x ,q x )g {q x ). (4.3) 

One should remind that when obtaining Eq. (|4.1I) we 
have omitted the average part of the kernel "E.{k x ,k' x ). 
The motivation to discard it, as was stated after 
Eq. H3.7|) . arises because its contribution to the renor- 
malization of the lengthwise wave number k n is real, i.e. 
it contributes only to the quantity j n (see Eq. I)l.l|l '). Be- 
sides this, another contribution to 7„ arises due to the 
real part of the self-energy. For this reason, we can ex- 
pand the average pole factor g(k x ) in series of partial 
fractions and retain only the imaginary part of the self- 
energy M(k n ), 



g(k x ) 



E 

n=l 



{■KU/df 



i 



i/2L f 



i/2L r 



(4.4) 



This expression completely corresponds to the represen- 
tation l|2.14c|) for the unperturbed pole factor. It is suit- 
able for further evaluation of the average Green's func- 
tion. 
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The quantity L n is the total wave attenuation length or 
electron mean-free-path of the n-th conducting mode. It 
describes the scattering from n-th mode into all possible 
propagating modes. This quantity is determined by the 
imaginary part of the self-energy M(k n ), 



k n d 



(4.5a) 



(Tm/d) 2 ^ (nn'/d) 2 



k n d k n > d 

n—l 

x[Q(k n ,-k n ,) + Q(k n ,+k n ,)}. (4.5b) 

Let us substitute Eqs. I|4.2a[) and (|4.4() into Eq. (|3.2() 
and perform straightforward calculations of the integrals 
with the use of the delta-function and over the poles of 
g(k x ). As a result, we find the average Green's function, 



(G(x, x'; z, z')) = G(\x - x'\;z, z'), (4.6) 



profile a^(x), and the latter from the terms depending 
on the roughness gradient o~£f{x). The expression (|5.2|) 
exactly coincides with that previously obtained by vari- 
ous methods (see, e.g., Ref. 

The second attenuation length L n ' related to RSGP 
spectrum through Q2(k x ,k' x ), is associated solely with 
the SGS mechanism due to the operator V(x) (see 
Eq. (|3.16a[n . In accordance with Eqs. (|3.24bl and 
(|3.21b|) . it is described by 



N d 



r(2) 



T (2) 

i'=l nn' 



(5.3) 



Here its diagonal term controlling the wave scattering 
inside the mode (intramode scattering) is written as 



L 



1 



a 4 (nn/d) 4 
2 kl 



1 



(2nn) 



x [T(2k n ) + T(0)}. 



(5.4) 



in the efficient representation via canonical Fourier series 
in the normal waveguide modes, 

Nd / ' \ 

g(\x - x'\;z, z') = sin i^j-) sin ( j 

n—l ^ ' 



exr>(ik n \x — x'l) / \x — x' 
x .. . exp 



i k"/i d 



2L 7 - 



(4.7) 



(2) 

The off-diagonal partial scattering length L±A n , that de- 
scribes the intermode scattering (from n-th mode to n'-th 
one, n ^ n'), is 

1 8a 4 (nn/d) 2 (irn'/d) 2 (n 2 +n' 2 ) 2 



r(2) 



7T 4 k„ 



k n , (n 2 - n' 2 ) 4 



[T(k n + k n ') + T(k n k n i 

)] ■ (5-5) 



ATTENUATION LENGTH ANALYSIS 



In view of Eq. (|3.25b|l , the general expression H4.5JI for 
the inverse attenuation length shows that in the problem 
under consideration this quantity consists of two terms, 



1 



-(i) 



(2)' 



(5-1) 



These terms descend from different mechanisms of the 
surface scattering. The first attenuation length Ln^ is 
related with the RHP spectrum, W(k x ), through the ex- 
pression for Qi(k x ,k' x ). In accordance with Eqs. Ij3.23bjl . 
(|3.18b|) and 13.19bjl . it is given by 

1 2 (im/d) 2 ^ (vrn'/d) 2 



-(1) 



k n d k n 'd 

n—l 



x [W(k n + k n >) + W(k n - k n >)] . 



(5.2) 



Its diagonal term is formed by the amplitude scattering 
(AS) and the off-diagonal terms result from the gradient 
scattering (GS). These two mechanisms of surface scat- 
tering are due to the corresponding terms in the expres- 
sion for S 1 (fc :E , k' x ) (see Eq. (JtOljaJ)), i.e., the former from 
the term depending on the amplitude of the roughness 



To the best of our knowledge, in the surface-scattering 
problem for multimode waveguides the operator V(x) was 
never taken into account, and, as a result, the second 

(2) 

attenuation length, or SGS length, L„ , was missed in 
previous studies. 

Now we list the simplifications that have been made 
in deriving Eqs. I|5.1|l - (|5.5|) for the attenuation length 
L n . First, the proper self-energy l|4.3|l in the Dyson- type 
equation for the average Green's function has been ob- 
tained within the second-order approximation in the per- 
turbation potential. In terms of the diagrammatic tech- 
nique this is similar to the "simple vortex" or, the same, 
Bourret approximation^ that contains the binary cor- 
relator Q(k x ,q x ) of the surface-scattering potential and 
the unperturbed pole factor ga(q x ). To find out the con- 
ditions of applicability for this approach, we have used 
the ideas proposed in the bookU More specifically, we 
substitute into the self-energy (|4.3|l the average pole fac- 
tor l|4.2b|) instead of the unperturbed one. This trick is 
equivalent to the summation of an infinite subsequence of 
diagrams in the exact expansion of the self-energy in pow- 
ers of the scattering potential. The analysis shows that 
in the Dyson-type equation the new (and more general) 
self-energy can be reduced to ours, if the channel broad- 
ening l/2L n is much less than the unperturbed quantum 
wave number k n and the variation scale R^ 1 of the RHP 
and RSGP spectra, i.e., when 2k n L n ^> 1 and R c <C 2L n . 
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Second, in order to extract the inverse attenuation 
length from the self-energy JOJ, we have changed the 
lengthwise wave number k x by its unperturbed value k„. 
This change is justified if the surface-induced broaden- 
ing l/2L n can be neglected in comparison with R~ l and 
the spacing \k n — k n ±\ \ « \dk n /dn\ between neighboring 
quantum wave numbers. Now we take into account that 
\dk n /dn\ ~ A" 1 , where A„ is the distance between two 
successive reflections of the n-th mode from the rough 
boundary. Therefore, the use of k n instead of k x in the 
argument of the self-energy is valid under the conditions 
R c < 2L n and A n < 2L n . 

Thus, we come to three requirements: R c -C 2L n , 
A„ <C 2L n , and 2k n L n ^> 1. Due to the obvious rela- 
tionship k n A n > 1, the last inequality is a direct conse- 
quence of the second one, and one can conclude that the 
domain of applicability for our results is restricted by two 
independent criteria of weak surface scattering, 



A, 



2k n d/(irn/d) < 2L Tl 



Rc < 2L r , 



(5.6) 



They imply that the wave is weakly attenuated on both 
the correlation length R c and the cycle length A„ . From 
the analysis performed above, it becomes clear that ex- 
pressions (|5.2|l and 1)5. 3[) - l|5.5ll represent main contri- 
butions from the substantially distinct surface-scattering 
mechanisms: AS+GS and SGS. In particular, the cor- 
rections that are proportional to ct 4 , originated from the 
next order of approximation in the amplitude and gra- 
dient terms of the surface-scattering potential, can not 
compete with the main contribution (j5.2|l under the con- 
ditions (|5.6|l . On the contrary, the square- gradient terms 
give rise to the cr 4 -terms in Eqs. I|5.4|l - l|5.5|l . which 
should not be neglected due to a specific dependence on 
the correlation length R c . Note that Eq. (j5.6|l implicitly 
includes the requirement a <C d for the surface corruga- 
tions be small in height, that has been employed in the 
section ITTT1 when deriving the explicit form (|3.6|l for the 
surface scattering potential S(fc x , k' x ). 

We are now in a position to analyze the attenuation 
length L n . For convenience, in our further analysis, we 
deal with the attenuation lengths in the form of the di- 
mensionless quantities A n /2L n , A n j2lff and A n l2lff ■ 
Since and iff depend on as many as four dimension- 
less parameters {ka) 2 , kR c , kd/ir, and n, the complete 
analysis appears to be quite complicated. For this rea- 
son, below we restrict ourselves by the analysis of the 
interplay between Ln and iff as a function of the di- 
mensionless correlation length kR c , for different values of 
{ka) 2 and mode index n. In the analysis, we consider a 
multimode waveguide, i.e., the situation when the num- 
ber of propagating modes is large, kd/n 3> 1. From 
the physical point of view, two types of rough surfaces 
seem to be the most important. Surfaces of the first 
type contain a small-scale roughness of the "white noise" 
kind when kR c <C 1. For the second type, the waveg- 
uide surface consists of large-scale random corrugations 
when kR c 3> 1. We shall develop our analysis for this 
two types of surfaces. 



A. Small-Scale Roughness 



Let us start with a relatively simple and widely 
used case of a small-scale boundary perturbation, when 
kR c <C 1 and the surface roughness can be regarded 
as a delta-correlated random process with the correlator 
W(x — x 1 ) W(0)5(x — x 1 ) and constant power spec- 
trum W(k x ) Ri W(0) ~ R c . Taking into account the 
evident relationship kA n > 1, one can get the following 
inequalities to specify this case 



kR c < 1 < kA n . 



(5.7) 



It is necessary to underline that in the regime of small- 
scale roughness l|5.7l) the second of the weak-scattering 
conditions in Ea. (|5.6[) is not so restrictive as the first 
one, and directly stems from it, R c <C A„ <C 2L n . 

In Eqs. (|5.2|) . (|5.4|l and (|5.5|l for the attenuation 
lengths the argument of the correlators W(k x ) and T(k x ) 
turns out to be much less than the scale of their decrease 
R^ 1 under the conditions (|5.7|l . Therefore, for any term 
in the summation over n' the argument can be taken as 
zero. 

Therefore, the first attenuation length is determined 
as follows, 



A TO 

2Ln ^ 



2(^7^™ Ma) 



kd/ir k k n >d 



(ka) 



o n kW(0) 



kd/n 



(5.8b) 



Due to a large number of the conducting modes ~ 
kd/n 3> 1, we can change the summation over n' by 



integration. In this way one can obtain Eq. I|5.8b(l from 
Eq. I|5.8a|) . In order to correctly estimate the result, one 
can take into account the formula 



dxW{x) = 2R C / dpW(Rcp), (5.9) 



W(0) 



which directly follows from the definition l|2.4|) for the 
Fourier transform W(k x ) of the binary correlator W(x). 
The function W(R c p) is the dimensionless correlator of 
the dimensionless variable p, with the scale of decrease of 
the order of one. As a result, the function W(R c p) does 
not depend on R c . Therefore, the integral over p entering 
Eq. I|5.9I) is a positive constant of the order of unity. For 
example, W(0) = \[2~nR c and the integral is \fnj2 in the 
case of Gaussian correlations (see Eq. I|5.27|) '). 
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TABLE I: Evaluation of Eq. lETTO 



n\An 


-4 


-3 


-2 


-1 


1 


2 


3 


4 


1 


0.002 


0.006 


0.024 


0.31 










2 




0.004 


0.019 


0.27 


0.31 








3 






0.018 


0.26 


0.27 


0.024 






4 








0.26 


0.26 


0.019 


0.006 




5 










0.26 


0.018 


0.004 


0.002 



For the SGS length we have 



A, 



; (fca) 4 



T(0) 



2L 



(2) 



k n d (fed/ 7r) A: 3 




Vn' — 1 n'— n+1 



7T 2 (k<j) 



16k n d 

■K^im/d) 2 

{nn'/d) 2 {n 2 + n' 2 ) 2 ' 
k n 'd 

T(0) 



(n 2 - n' 2 ) 4 



4 fe n d (kd/Tr) k 3 



(5.10a) 



(5.10b) 



In Eq. I|5.10af> every term in the sum rapidly decreases 
with an increase of absolute value of An = n — n' . This 
can be seen by making use the following estimate, 

(n 2 



a\2 



l 



(n 2 - n' 2 Y 4(An) 4 



for n ^> |An| 



(5.11) 



The fast decrease of the factor is supported by direct cal- 
culations, see Tabled Therefore, the sum in Eq. I|5.10a|) 
can be well evaluated by two terms with n' = n ± 1. For 
simplicity, in Eq. ((5.1 Obi) we assume AT d ^> n 3> 1, and 
replace the curly braces by factor 1/4. 

The explicit form for T(0) directly follows from the 
definition 1(3.13(1 for the correlator T(k x ), 



/oo 
dxW" 2 (x) 
-oo 



= R r 



dp 



d 2 W(R cP ) 
dp 2 



(5.12) 



If the roughness correlations are of the Gaussian form, 
then according to Eq. (JS2HI, we have T(0) = 3y/n/iRl 
and the integral over p entering Eq. ((5.121) is equal to 
3yW4. 

The relationship between the attenuation lengths is ex- 
pressed by 



r(i) „2 



k n d 



(5.13) 



According to substantially different behavior of the quan- 
tities A n /2Ln^ and A„/2L^ 2 ' > with respect to kR c , it be- 
comes clear that they must intersect at the crossing point 
(kR c ) cr . If the crossing point falls onto the present re- 
gion of small-scale roughness (kR c <C 1), its dependence 
on the model parameters is obtained by equating to one 
the expression 1(5.13(1 . 



(5.14) 



To the left from this point (kR c ) cr the SGS length pre- 
vails, Ll 2) < 1$. To its right the main contribution 



is due to the first attenuation length, L n 



(i) 



The 



expression (15.14(1 shows that the crossing point is smaller 
for smaller values of the dimensionless roughness height 
ka, as well as for smaller mode indices n, or for larger 
values of the mode parameter kd/ir. 



B. Large-Scale Roughness: Weak Correlations 

The intermediate situation arises when the correlation 
length R c becomes much larger than the wave length 
2it/k, but still remains much less than the cycle length 
A„, 



1 < kR c < fcA, 



(5.15) 



As before, the first of the weak-scattering conditions 15. 6|) 
is the most restrictive, i.e., we get R c <C A„ <C 2L n . 

Since the distance A„ is larger than the correlation 
length R c , successive reflections of the waves from the 
rough surface are weakly correlated. Meantime, the dis- 
tance between neighboring wave numbers k n and k n ±\ is 
much smaller than the variation scale R^ 1 of the corre- 
lators W(k x ) and T(k x ), 



\k n - k n±1 \ » \dk n /dn\ = 27TA- 1 < R-\ 



(5.16) 



This implies that the correlators W(k x ) and T{k x ) are 
smooth functions of the summation index n' . Therefore, 
the sum in the expression 1(5.2(1 for the first attenuation 
length Ln can be substituted by the integral, 



A, 



2L 



(i) 



dn 



,{™>/d) 2 



x [W(k n + k n ,) + W{k n -k n ,)] 



a /7rn\ 
= ~ \~d) 



dk x ^k 2 -k 2 x W{k n -k x ). 



(5.17a) 
(5.17b) 



Eq. 1(5.17(1 shows that the first attenuation length is con- 
tributed by scattering of a given propagating n-th mode 
into all other propagating modes. Note that to obtain 
this asymptotic result we have used only the condition 
of week correlations, R c <C A n . Therefore, Eq. ((5.17|) 
provides the reduction to Eq. 1(5.8(1 for small-scale corru- 
gations, kR c <§; 1. In the case of large-scale roughness, 
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kR c 1, the formula (|5.17() obtained for allows 
further simplifications as was done in Ref. U 

In contrast to L n x \ due to a rapidly decaying factor 
(|5.11|) . the SGS length L%> can be still described by 
keeping tree terms only, n' = n, n ± 1, in the sum in 
Eq. (|5.3|) . Taking into account the estimate 15.16|l . for 
the case Na ^> n^> 1 one can write down, 



A„ 


_ 7T 2 (fcfj) 4 


n 3 


T(0) + T(2k n ) 




2L n 


) '~ 2 fc„<2 


(kd/n) 


k 3 




■I 


'1 1 


2 8 


} 


(5.18a) 


3 + (27m) 2 _ 






_ tt 2 (fed) 4 n 3 


r(o)- 


h r(2fe„) 


(5.18b) 


8 


k n d (kd/n 


) 


fc 3 



In the final expression (l5~18bl) . we have replaced the curly 
braces from Eq. (|5.18a|l by the factor 1/4. Naturally, at 
small-scale corrugations, kR c -C 1, the obtained result 
(|5.18l) passes into Eq. (|5.10|) . For the large-scale rough- 
ness, when kR c ^> 1, one should use Eq. (|5.18() because 
of arbitrary value of the parameter k n R c . We do not 
consider here this case in detail due to its intermediate 
character. 



C. Large-Scale Roughness: Strong Correlations 

In the other extreme case, the correlation length R c is 
very large not only in comparison with the wave length 
27r/fc, but also in comparison with the cycle length A„, 



1 < kA n < kR c . 



(5.19) 



In this case the number of wave reflections over the cor- 
relation length R c is large. Therefore, the successive re- 
flections are strongly correlated to each other. 

Under the inequalities (|5.19|) the second of the weak- 
scattering conditions l|5.6|) is the most restrictive, there- 
fore, the condition of applicability reads as 



A„ < R c < 2L r . 



(5.20) 



The latter requirement determines the upper limit for the 
value of the correlation length R c . 

Due to Eq. (|5.19|) . the distance between neighboring 
wave numbers k n and k n ±\ turns out to be much larger 
than the variation scale R~ l of the correlators W(k x ) 
and T(k x ), 



\k n fcn±l| 



\dk„/dn\ = 27TA- 1 > R~ x . 



(5.21) 



This indicates that the probability of the intermode 
(n' 7^ n) transitions is exponentially small and the atten- 
uation lengths are mainly formed by the incoherent in- 
tramode (n' — n) scattering. Formally, at strong correla- 
tions l|5.19|l the correlators W{k x ) and T(k x ) are sharpest 
functions of the summation index n' . In the sums of 



Eqs. I|5.2(l and l|5.3|) for the attenuation lengths the main 
contribution is due to the diagonal terms with n' = n, in 
which W{2k n ) and T(2k n ) can be neglected in compari- 
son with W{0) and T(0). 

Thus, the first attenuation length reads 



A„ (kaf 



2Ln ^ 



kW{0). 



k n d (kd/n) 3 
Correspondingly, for the SGS length one gets, 
7T 2 (fcfi) 4 n 3 T(0) 



A„ 
2il 2) 



2 k n d (kd/n) k 3 



1 



1 2 



(27m) 2 



(5.22) 



(5.23) 



The ratio of the first attenuation length to the second 
one can be presented as 



r(l) 

7W) 



A, 



2L 



(i) 



(fc„A n ) 2 < 1. 



(5.24) 



According to the first inequality in Eq. (|5.6(l , to the con- 
dition (|5.19|l of the strong-correlations and to the evi- 
dent relationship fc„A„ > 1, we can see that the am- 
plitude scattering length always prevails over the SGS 
length within the interval of strong correlations. For this 
reason, in the condition of applicability H5.20JI one should 
substitute L n by Ln ■ This allows one to arrive at the 
inequalities in the explicit form, 



kA n < kR c < (kA n )(kcr) 1 (kd/im). 



D. Numerical Analysis 



(5.25) 



In this subsection we perform the numerical analysis 
of the scattering lengths for the case when the random 
surface profile £(x) has the Gaussian binary correlator, 

W(x) = cxp(-a; 2 /2 J R 2 ). (5.26) 

Then the RHP spectrum (|2.4() is given by 

W{k x ) = V2ttR c exp(-fc 2 i? 2 /2). (5.27) 

The RSGP spectrum T(k x ) defined by Eq. (|3~T3|) ). can 

be explicitly presented as 



T(k x ) = 



[{k x R c f - A{k x R c ) 2 + 12] 



l&Rl 
exp [-(k x R c ) 2 /4] 



(5.28) 



In Fig. ^ we display separately the behavior of 
A n /2Ln and A n /2Ln \ as a function of the dimension- 
less correlation parameter kR c , comparing them with the 
analytically obtained asymptotics. The dashed lines are 
used to plot the asymptotics (|5.8|l and (|5.10|) for the re- 
gion H5.7fl of small-scale roughness (kR c <C 1 < kA n ), 
and expressions l|5.22|l and l|5.23|l for the region (|5.19|) 
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kR 

c 

FIG. 1: Plots of A n /2L^ (increasing curves) and A n /2L^ 
(decreasing curves) vs. kR c for kd/n — 100.5. Dashed lines 
show the corresponding asymptotic expressions. 



of large-scale roughness with strong correlations (1 < 
kA n -c kRc). As one can see, for these two regions nu- 
merical data for both lengths are quite well described by 
the corresponding asymptotic expressions. 

The curves A n /2£„ clearly manifest two transition 
points: between the regions of small-scale and large-scale 
corrugations at kR c ~ 1, and between weak and strong 
correlations at kR c ~ kA n . As follows from Eq. I|5.2jl . 
the inverse value of the first attenuation length typically 
increases with an increase of kR c . Specifically, within the 
interval of the small-scale roughness (kR c <C 1 < kA n ) we 
have 1/Ln ^ kR c . Then, within the intermediate region 
of large-scale roughness with weak correlations where 
1 <C kR c <C fcA„, the increase of A n /2L„' slows down 
(see the curves with the parameters (ka) 2 = 10~ 2 , 10~ 4 ; 
n = 50), or can even be replaced by the decrease for some 
values of the model parameters (see the curve with the 
parameter (ka) 2 = 1CP 2 , n = 10). Within these two re- 
gions that are unified under the condition kR c <C kA n , 

the quantity is determined by both AS and GS 

mechanisms (AS+GS). Finally, for large-scale roughness 
and strong correlations (1 < kA n <§; kR c ) the value of 

again begins to increase linearly with kR c . Here 

is associated solely with the AS mechanism be- 
cause the main contribution to the asymptotic i|5.22[l is 
due to the diagonal term in the sum i|5.2|) . 

In contrast with 1/ Ln \ the inverse SGS length l/lffl 
reveals a monotonous decrease as the parameter kR c in- 
creases. At small (kR c <Cl< kA n ) and extremely large 
(1 < kA n <C kR c ) values of kR c , this decrease obeys the 
law 1/Ll 2) cx {kR c )- 3 , due to T(0) - R~ 3 . 

In Fig. n we can also see the crossover from the SGS 
to AS+GS that is characterized by the crossing point 



III 1 

Jimall-scale- [ Large-scale-roughness 

■ \ \V_ ^ / 
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c 

FIG. 2: Plots of A n /2L„ versus kR c for fcd/vr = 100.5, 
n = 50, and different values of (ka) 2 . The dot-dashed lines 
labelled '(fc_R c ) cr ' and 'kAn, indicate the transition between 
the regions dominated by the SGS, AS+GS and AS mecha- 
nisms (from the left to the right). The dashed line at kR c = 1 
marks the transition between the regions of small-scale- and 
large-scale-roughness . 



between A n /2L ( n > and A n /2L ( n >. The crossing point for 
two curves with the parameters (ka) 2 = 10~ 4 and n — 50 
is very close to that for two curves with the parameters 
(ka) 2 = 10~ 2 and n = 10. Approximately, both crossing 
points are 

(kR c ) cr - 0.2. (5.29) 

They are located well inside the interval of small-scale 
roughness, and their values are in an agreement with the 
asymptotic expression (|5.14|l . Two curves corresponding 
to the parameters (ka) 2 — 10~ 2 and n — 50, have cross- 
ing point in the transition region kR c ~ 1 between small- 
and large-scale corrugations. 

In the following two Figures we display the dependence 
of A n /2L n as a function of kR c . The curves are plotted 
starting from the values of kR c for which A„/2L„ — 1, 
according to the first condition in l|5.6|l . Taking into 
account the second condition restricting the maximal 
value of kR c , we plot every curve within the range where 
R c < 2L ( n } . Based upon the description of the Fig.^ the 
identification of each scattering mechanism dominating 
in the corresponding regions becomes simple. One can 
see that the curves in Figs.[21and[!5]cxperience firstly the 
crossover from the SGS to the AS+GS, and after, from 
the AS+GS to AS. We outline both transitions with the 
labels l (kR c ) cr ' and 'fcA„', respectively. 

In Fig. |21 we show the dependence of A n /2L n on kR c 
for three values of the parameter (ka) 2 (two of them 
(ka) 2 = 10~ 4 , 10 -2 correspond to those used in Fig.^l. 
The curve with (ka) 2 = 10 -4 has the crossing point 
(kR c ) cr of the value l|5.29|l located within the interval 
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FIG. 3: Plots of A n /2L n versus kR c for kd/iv = 100.5, 
(ka) 2 = 10 -2 , and different values of the mode index n. The 
sets of square symbols labelled '(fci? c )cr' and 'fcA„' play the 
same role as dot-dashed lines in Fig. [5] 



of small-scale roughness. The crossover reveals a small 
dip centered at {kR c ) cr . The curve obeys the asymptotic 
behavior (fci? c )~ 3 to the left from (kR c ) cr due to the 

(2) 

main contribution from A n /2L n . After, the quantity 

A n /2Ln^ becomes dominating in the sum l|5.1|) . there- 
fore, the curve begins to rise. Firstly, the linear depen- 
dence on kR c on the right deep-side (where kR c < 1) 
is replaced with a smoother one (for kR c > 1). Finally, 
for R c > A n (strong correlations) the linear dependence 
restores. 

The crossing points of the second and third curves 
with (ka) 2 = 10 _2 ,10 _1 have values of the order of 
unit, (kR c ) cr ~ 1. Here the total attenuation length L n 
within the whole small-scale region is formed solely by the 
SGS length L { n\ In full agreement with Eq. (|5~l"4) the 
presented curves display that the smaller the parameter 
(ka) 2 the smaller the value of the crossing point (kR c ) cr . 

To visualize the dependence on the mode index n, in 
Fig. |2| we plot A n /2L n with the parameter n spanning 
the range of the propagating modes. One can see that for 
the mode index n = 1 the crossing point (kR c ) cr ~ 0.1, 
the curve with n = 10 has (kR c ) cr given by Eq. 15.29fl . 
For the rest of the values of n the crossover occur at 
(kR c ) cr ~ 1.0. We note that the smaller the channel in- 
dex n, the smaller the value of the crossing point (kR c ) cr . 
This fact is in full agreement with Eq. I|5.14|l . The squares 
labelled with l kA n \ show the transition points between 
the AS+GS and AS mechanisms (or between the weak- 
and strong-correlations of the surface roughness). They 
are fcAi » 63500, fcAi » 6315, fcA 50 m 1100, fcA 90 « 310 
and fcAioo ~ 60. 

Finally, note that for all curves in our figures the 
roughness height is small, a/d <C 1. Furthermore, for 
the amplitude- and gradient-dominated scattering (to the 



right from the point (kR c ) cr where A n /2Ln mainly con- 
tributes), the average corrugation slope is also small for 
all data, a/R c -C 1. The roughness slope remains to be 
small at the crossing points too, but increases to their left 
with the decrease of kR c . As a result, to the left from the 

(2) 

crossing point where the square-gradient term A n j2L n 
prevails, the slope reaches the values of the order one, or 
even larger. 



VI. CONCLUSION 

In this paper we investigated the wave/electron scat- 
tering in multi-mode surface-corrugated waveguides of 
quasi- ID geometry. For this kind of waveguides, we have 
discovered a new square- gradient scattering (SGS) mech- 
anism that is different from the previously studied ones. 
This mechanism arises due to specific square-gradient 
terms in the Hamiltonian describing the surface scat- 
tering, that are related to the roughness square-gradient 
power (RSGP) spectrum T(k x ). To compare with, the 
well known scattering mechanisms, the amplitude (AS) 
and the gradient (GS) ones, are both determined by the 
roughness-height power (RHP) spectrum W(k x ), only. 
Since the SGS mechanism is independent from the oth- 
ers, one can define two attenuation lengths, the known 
length Ln and the SGS length Ln \ Both contribute to 
the total attenuation length (or, the same, electron mean 
free path) L n according to Eq. jOJ . 

The roughness-height W(k x ) and square-gradient 
T(k x ) power spectra have very different dependencies 
on the roughness correlation length R c . This provides 
the substantially different behavior of the corresponding 
scattering lengths in dependence of the model parame- 
ters. Specifically, the inverse value of the first attenua- 
tion length 1/Ln typically increases, while the inverse 

(21 

value of the SGS length 1/L n decreases, with an in- 
crease of the parameter kR c . Therefore, the curves dis- 
playing these quantities intersect upon the increase of the 
dimensionless correlation length kR c , and the crossover 
from the SGS to AS+GS occurs. To the left from the 
crossing point (kR c ) cr the SGS length prevails over the 

first attenuation length, Ln <C L„- To the right from 
(kR c ) cr , the first attenuation length mainly contributes 
to the scattering process, Ln *C Ln ■ If the crossing 
point (kR c ) cr falls into the interval (|5.7|l of small-scale 
surface corrugations, it obeys the law (|5.14|l . 

As we have shown, at any fixed value of the root- 
mean-square roughness height a, one can indicate the 
region of small values of the correlation length R c where 

(2) 

the new attenuation length L n predominates over the 
known length L„ . This predominance arises in spite of 
the fact that l/L n is proportional to a 2 while 1/L^ is 
proportional to c 4 . 

In the large-scale roughness regime where the first at- 
tenuation length mainly contributes, L^ <C , one 
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can observe two different behaviors of \/Ln\ In the tential 
interval of weak correlations (|5.15|) the dependence of 
on kR c is quite complicated, due to the coex- 
istence of both the AS and GS mechanisms. How- 
ever in the region of strong correlations (|5.19() . because 
the AS stands alone, we have much simpler behavior, 



oc kR c . It is remarkable that the SGS mecha- 
nism prevails in the widely discussed region of a small- 
scale boundary perturbation, kR c <C 1, where the surface 
roughness is typically described via the white-noise po- 
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